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DIRECT METHOD OF DESIGN AND STRESS ANALYSIS OF ROTATING 
DISKS WITH TEMPERATURE GRADIENT 

By S. S. Masson 


SUMMARY 

A method is •presented for the determination of the contour of 
disks, typified by those of aircraft gas turbines, to incorporate 
arbitrary elastic-stress distributions resulting from, either centrif- 
ugal or combined centrifugal and thermal effects. The specified 
stresses may be radial, tangential, or any combination of the two. 
For example, the equivalent stress may be specified at each radius 
inasmuch as this stress is generally used as the criterion of the 
occurrence of plastic flow. Use is made of the finite-difference 
approach in solving the stress equations, the amount of computa- 
tion necessary in the evolution of a design being greatly reduced 
by the judicious selection of point stations and by the aid of a 
design chart. Use of the charts and of a preselected schedule of 
point stations is also applied to the direct problem of finding the 
elastic- and plastic-stress distribution in disks of a given design, 
thereby effecting a great reduction in the amount of calculation 
compared with previously published methods. Illustrative 
examples are presented to show computational procedures in the 
determination of a new design and in analyzing an existing 
design for elastic stress and for stresses resulting from plastic 
flow. 

INTRODUCTION 

The important role of the disk in the gas turbine has 
directed much attention to the analysis of stresses in rotating 
disks with temperature gradient. Early in the development- 
of stress analysis of steam turbines, Stodola (reference 1) 
presented the basic equation for stresses in such disks, as 
well as a limited approach toward the solution. More 
recently, in connection with the development of the gas 
turbine for aircraft propulsion, the problem of determining 
thermal stresses resulting from temperature gradient has 
been actively investigated. Several methods have been 
developed (references 2 to 5) that can be used to calculate 
the stresses in a given disk due to differential expan- 
sion and to variation in physical properties of various 
parts of the disk operating at different temperatures. Thus, 
for disks already designed, it is readily possible to determine 
with a high degree of accuracy the centrifugal and thermal 
stresses resulting from a- given speed and temperature 
distribution. 

Little attention has, however, been directed towards the 
problem of designing the profile of a disk, which, when oper- 
ated at a given speed and temperature distribution, will 
incorporate a given total-stress system— centrifugal plus 


thermal. A direct method of design for prescribed total- 
stress systems was developed at the NACA Lewis laboratory 
during 1948. Although the design method can readily be 
extended to include plastic flow, in the interest of simplicity 
only elastic stresses are considered. 

The specified stresses may be either radial or tangential 
values, as is common in engineering practice, or the speci- 
fication may refer to any combination of these stresses. For 
example, specification of the equivalent stress at each radius 
may be desired inasmuch as this stress is frequently used as 
a criterion of the occurrence of plastic flow. The method is 
based on the finite-difference equations derived in reference 3 
and makes use of several simple design charts, presented 
herein, that facilitate the calculations. An illustration is 
presented for tracing the actual steps that may be followed 
in evolving a design. 

A simplified method of analyzing an existing disk for 
centrifugal and thermal stress is also described. The method 
is especially useful in checking the stresses in a disk designed 
for prescribed stress distributions, particularly if min or 
deviations from the design profile are introduced for practical 
reasons, or if creep or plastic flow, not considered during the 
design stage, are likely to occur. Applications to arbitrary 
disks are equally valid. Essentially, the method is an 
adaptation of the finite-difference approach of references 3 
and 4, the principal modifications lying in the use of a pre- 
selected distribution of point stations and the addition of 
several alinement charts to facilitate the computations. 

Although the methods presented are of special advantage 
in calculations involving thermal stresses, they are suitable 
for disks subjected solely to centrifugal effects. Simple 
equations may be obtained both for design or for stress 
analysis of rotating disks by setting all terms involving tem- 
perature equal to zero. The subject of design for prescribed 
centrifugal-stress distribution has likewise received only 
limited attention (examples of which may be found in refer- 
ences 6 and 7) ; the present method may therefore satisfy a 
need in this field of design. 

DIRECT METHOD OF DISK DESIGN 

Equations. — The equations for stress analysis of rotating 
disks are (reference 3) 

^ (rAo>) — htr t -r p(Ahr s =Q (1) 
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and 

d {<T t \ d /po>\ , d (1+m)(o>— cr,)_ n /nv 

dr \e) dr \E ) + dr {aAT) ' Er “° (2) 

(All symbols are defined in .tie. appendix.) 

Equation (2) is a statement, of the relation that must exist 
between the radial and tangential stress distributions, 
regardless of the shape of the disk, and therefore does.. not 
conjtain the disk thickness h. Once the stresses have been 
chosen to be compatible in accordance: with equation (2), h 
can be determined from equation (1). 

Solution of compatibility equation for specified stresses. — 
A finite-difference method for the solution of equations (1) 
and (2) is presented in reference 3. Choosing discrete point 
stations along a radius permits the stresses at the._«th 
station to be expressed in terms of the stresses at the (n-l)st 
station. Equation (2) in finite-difference form becomes 

C'n<r r ,n— — G'nfft.n-I+H'* (3) 

The coefficients C r n , DC, and so forth are, in general, com- 
plicated expressions involving the radii r n and r„_ i at the 
nth and (n-l)st stations as well as the properties of the 
disk material at the two stations. A great simplification 
in the expressions results if the radii r„ and r B _i are not 
chosen at random, but are chosen according to a preselected 
schedule, for example, according to the relation 

0,200 

r n 

If, in addition, Poisson’s ratio p is assumed constant and 
equal to 0.333, equation (3) reduces to 

1.133 <r liB — 0.466 v Tin =-M^~ (0.833 o-(, B _i— 0.166 o>,„_i) —EJEI' n 

Mhn-i 

=Z* (4) 

Equation (4) is useful in the central region of the disk 
where values of r K are relatively low. In order to obtain 
reasonable spacing of the point stations in the rim region 
where the radii are large, it is desirable to choose the spacing 
according to the relation 

r *— r "- 1 =Q .050 

7 1 n 

Equation (3) then becomes . . _ 

w 

1.033 <r*, B — 0.366<r r , n =-^r^- (0.965ff, B _i — 0.298ffr.»-i) — 

-C'n-l 


specified^ In some designs, specification of either the radial 
or the tangential stress may be desired; hence, the unspecified 
stress may then be directly determined from equation (4) or 
from the corresponding straight-line plot. A disadvantage of 
starting a design by specifying the radial or tangential stress 
is that the unspecified stress may be too high or in some 
other way unfavorable. It may be more desirable to specify 
the equivalent stress <r e , which is defined by the equation 

O’e.n — V°V,n* l,ii* (0) 

This stress is usually used as the criterion for plastic-flow 
calculations (reference 4) and is therefore a logical design 
parameter. Furthermore, specification of the equivalent- 
stress insures that neither the radial nor the tangential stress 
will be unexpectedly high. A plot of <r Till against er t> „ for a 
given value of <r«,„ results in an ellipse inclined at 45° to the 
coordinate axes. The intersection of the line defined by 
equation (4) and the ellipse of equation (6) defines the specific 
values of o> ( „ and <x t-n . Although, in general, a straight line 
intersects an ellipse at two points, the proper point for 
design use is usually obvious from the physical nature of tho 
problem. 

Design chart. — A design chart based on equations (4) and 
(6) is presented in figure 1(a). The straight lines correspond 
to different values of Z a in equation (4) and the ellipses 
correspond to different values of equivalent stress This 
chart is used to calculate the stresses at the nth station after 
the stresses have been determined at the (n-l)st station. 
The value of Z n is calculated by equation (4) from conditions 
at the (n-l)st station and the temperature at the nth station, 
thereby establishing the line along which o>, B and <r t , K must 
he. Intersection of this lino with the ellipse corresponding 
to the assigned value of equivalent stress at the nth station 
establishes cr r;B and cr ( ,„. Thus, starting at a station near tho 
center of the disk (or for a disk with a central hole, at tho 
edge of the hole), where the radial and tangential stresses 
bear known relations to the equivalent stress, makes it 
possible to proceed with the determination of stresses at 
successive stations of increasing radius. 

The corresponding chart for equation (5) is figure 1(b), 
which is used in determining the stress conditions in tho rim 
region; whereas figure 1(a) is used for determining the stress 
conditions in the central region of the disk. 

Determination of disk profile. — When the stresses have 
been determined, the shape of tho disk can be obtained from 
the finite-difference equivalent of equation (1). For random 
values of r n and r n _ 1; the equation assumes the form 


E„H\—Z n (5) 

As an example, equation (4) is considered. If <r r , B _rand 
<r ( , B _ i are already established and if the temperatures at the 
nth and (n-l)st station are known, the value of Z n can be 
determined; then the only unknowns are <r r , B and o- t>B . The 
equation is linear in <r T , n and <r ( ,„ so that whatever are the 
actual values, a plot of o>,„ against v {>B results in a straight 
line defined by equation (4). In order to establish the 
individual values of <r T , n and c !iB , another condition must be 


h ( r -r ) ff r,n-Vi,*+p s wV 

tT / 2r„_! \ I - ' 

Z~~ ) O'r.s-l+O'f.B-l — 

V n r n~ls 


For — 


= 0.200 and pn—pn-i = constant density, equa- 


tion (7) becomes 


h n - 1 _ 10.00g r , B — g- ( , B +p B r B y 

h a 8.00o>, B _i+ff< r »-i— 0.640p„r B s c<) 4 


( 8 ) 



Tangential stress, cr, )n , Ib/sq in. 
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and for — =0.050 and p fl =p„_i=constanfc density, 

r n .... .... .... . 

equation (7) becomes 


presumably has already been determined from the disk 
design. Hence, 

Total resisting force at rim of continuous disk 


40.00gy,tt ^ PrJ'n^ /g\ 

h n 38.00o - r,*-i _ l _ & t,n-\ 0.903p s r a 8 w 3 

Because equations (8) and (9) give only ratios of thick- 
ness and not specific thicknesses, a starting point for calcu- 
lating actual thicknesses must be obtained from other 
considerations, for example, those at the rim. For disks 
with the familiar fir-tree blade attachments, one approach 
is as follows: 

With reference to figure 2, if sufficient clearance is provided 
in the fastening between the blades and the disk to prevent 
circumferential tightening, the radial stress in the disk at the 
base of the serrations (at radius r b ) is the sum of the cen- 
trifugal stresses produced by the blades and by the serrated 
regions of the rim. 



Figure 2. — Disk nomenclature for derivation of rim thickness. 




The total centrifugal force due to airfoil sections of blades is 

, NffbiAti ( 10 ) 

The serrated bases of the blades and the base-retaining 
sections of the disk may be considered as a discontinuous ring 
of material subjected to centrifugal effect, the serrations 
serving to eliminate tangential stresses. When it is assumed 
that at the rim the thickness h d (fig. 2) is specified from con- 
siderations of blade design and that the thickness tapers 
linearly to a value h b (as yet to be determined) at the base of 
the serrated section of the disk, then approximately: 

Total centrifugal force on .serrated sections of blades and 
disk is 

(speed 2 ) (volume) (average density) (average radius) 
or - 

»-(r ( -r0 2„ p ^ (11) 


= (stress) (area) 

■ =ff Tit ,2Tri,hi > ( 12 ) 

When the sum of equations (10) and (11) is equated to 
equation. (12), 

, 47Sfo- M ^ i , i +vpM 2 (r g +r t ) t (r d — r b )h d 

* Sttr b a T , b — 7rpw s (r s +r ll ) 2 (r tf — r s ) 


When h t has been determined, the thicknesses at the remain- 
ing stations in the disk may be obtained by successive 

multiplication by the ratio An application of the 

method is illustrated by the following example. 

ILLUSTRATIVE DESIGN 


Requirements. — -The assumed problem is to design a disk 
of 9-ine.h radius to carry blades that require a 1.9-inch base 
thickness (h d = 1.9). In order to retain the blades, a serrated 
region of 1-inch radial depth is necessary (r t = 9-1 = 8 in.). 
Based on measurements of similar disks, the temperature 
distribution is assumed to be as shown in figure 3, and as an 
initial design, the equivalent stress at each location is assumed 
to be equal to the elastic limit of the material at the operating 
temperature. The disk is to operate at a speed of 1 1 ,850 rpm 
and the density of the disk and blade material is 0.29 pound 
per cubic inch. 



Radial distance, r, in. 

Figure 3.— Variation of temperature with disk radius assumed for Illustrative design. 


Establishment of stresses satisfying the compatibility 
equation. — The calculations are shown in table I. The disk 
is divided into 18 stations. Column 1 indicates the radius at 
each station as a fraction of the radius of the continuous 
section of the disk. These values have been chosen to yield 

nine stations with - ■■ '*" ■=0.200 and nine stations with 



=0.050. 


The actual value of r n at each station is 


The radial stress at the base of the slots is equal to the radial 
stress ov.6 at the rim of the continuous section of the disk and 


r, 
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obtained by multiplying the values of column 1 by the disk 
radius, which in this case is 8 inches. The radii are listed in 
column 2. 

The difference between the actual temperature and the 
temperature at which there is no thermal stress is listed in 
column 3. Because in this case the stress-free condition is at 
a room temperature of 70° F, the A T n term of column 3 is 
obtained by subtracting 70° F from the actual temperatures 
in figure 3. 

The elastic moduli and the coefficients of thermal expan- 
sion at the operating temperatures of the various stations are 
listed in columns 4 and 5, respectively. 

The values a n AT n , and E„H* 

are then calculated as shown in columns 6, 7, and 8, respec- 
tively. 

The quantity r x 2 p n a? at each station is then calculated, as 
shown in column 9. Column 10 lists values of the ratio 


E n 


needed for subsequent calculations. 


Columns 11 and 12 are temporarily bypassed, and in col- 
umn 13 the desired equivalent stress at each station is 
entered. In this design, the equivalent stresses are taken 
equal to the elastic limits of the material at the operating 
temperatures, which for the temperatures of figure 3 are 
shown in column 13. 

The essential stress determinations are then effected by 
simultaneous operation on columns 11 to 15. At station a, 
the radial and tangential stresses are both equal to <r, and are 
therefore both entered as 76,500 pounds per square inch in 
columns 14 and 15. These stresses are used to determine 
the value of Z at station 2 equal to 49,600, as listed in column 
12. Reference to figure 1(a) then shows that the line Z= 
49,600 intersects the ellipse <r,=76,300 at o- r ,„=76,000 and 
at cr tiW =75,000, which are then listed in columns 14 and 15. 
In the same maimer, these values are used to determine the 
stresses at station 3 and progressively at subsequent stations. 
For example, at station 8, Z=32,200 and <r«=75,000. Inter- 
section of the corresponding straight line and ellipse yields, 
in figure 1(a), <jv.*=83,000 and <r,,*=62,000. 


At station 10, there is a change to — — =0.050. The 

7 * * 

only effect is to change the coefficients in column 1 1 and to 
introduce the use of figure 1(b) instead of 1(a). At station 
10, for example, Z— 13,800 and o-,= 73,000. Intersection of 
the corresponding straight line and ellipse in figure 1(b) 
yields o>,*=83,500 and 42,500, which are listed in 
columns 14 and 15, and are subsequently used to obtain 
the value for Z at station 11. Thus, by successive and 
simultaneous operation on the various columns, values of 
<r Ti% and may be obtained for each station that are con- 

sistent with the compatibility equation and combine to yield 


an equivalent stress at each station of any desired value. 
(If the imposed requirements are impossible to attain, the 
necessary straight line and ellipse will not intersect. Thus, 
for example, were it arbitrarily required to achieve an equiv- 
alent stress at the rim of 40,000 lb/sq in. while maintaining 
the equivalent stresses at all other stations at the values listed 
in column 13, the line Z= 60,900 would not intersect the 
ellipse tr,=40,000. The indication is merely that no physical 
disk can satisfy the requirements and that the requirements 
must be changed.) 

Determination of disk profile. — Establishment of stresses 
by means of the compatibility equation makes it possible to 
obtain thickness ratios in accordance with equations (7) and 
(8). Before proceeding, however, examination of the radial 
stress at the rim ay t!) should be made in order to determine 
whether or not it will yield a reasonable value of thickness h b 
in accordance with equation (13). Substitution of the values 

iV=54 (rj{=27,000 

-A 4 ,=0.625 p&!*=1150 

r s =8.00 r d =9.00 


A*=1.9 


in equation (13) results in the relation between h„ and a- rfi 
shown in figure 4. Obviously, from this figure, values of 
<r r .j below 18,000 pounds per square inch will yield inordi- 
nately large thickness values of k b , whereas values of <r rJ> 
above 40,000 pounds per square inch will yield inordinately 
low values while overstressing the rim material, which is 
operating at a high temperature. The allowable radial rim 
stress can be ascertained from a detailed analysis of creep 
strength and stress-rupture values of the rim material for 
the desired life of the wheel; but in any case, figure 4 indi- 
cates that the reasonable range for is between 18,000 
and 40,000 pounds per square inch. 



Figure 4.— Relation between thickness at rim of continuous section of disk and radial stress 

at that location. 
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Obviously, the value of a T , h in table I obtained in the first 
determination (o>,&= approximately 0) does not lie within 
the aforementioned reasonable limits. The reason for the 
very low value of <r r>s , however, is the designation of the low 
value of o’,, 6 =58,000 pounds per square inch at rim of con- 
tinuous section of disk. The compatibility condition 
demands only that <r t ,b and <r t , t lie along the line Z= 60,900 
(column 12, table I). By modification of the demand .on 
cr e ,6, the individual values oi c t .b and cr r , 6 can be moved to 
other locations along Z=60,900. Thus, for example, by a 
change in the requirements on a e , b to 63,500, 64,000, 67,700, 
and 72,900 pounds per square inch, the values of <j t , » can be 
successively moved to 18,00.0, 20,000, 30,000, and 40,000 
pounds per square inch with corresponding changes in <r t y, 
all these values are given in table I. Physically, the signifi- 
cance of this result is the impracticability of attempting to 
find a disk for which stresses, At the rim will not exceed the 
elastic limit if all the other, regions are to operate exactly at 
the elastic limit. Allowing the rim to operate somewhat 
above the elastic limit permits the remainder, of the disk to 
be kept at the elastic limit. Under, unusual conditions, 
modifying several stations near the rim may be necessary. 


The values of 


are now computed according to equa- 


tions (8) and (9) , the necessary calculations being shown in 
columns 16, 17, and 18. Actual values of h n are obtained 


in column 19 by successive, multiplication by the ratio 


h n -i 

K 


starting with the rim value h b determined from figure 4 for 
each value of v T%b and proceeding toward the center of the 
disk. The disk profiles and the stress distributions are 
shown in figure 5. The slenderness of the profiles relative 
to the rim is the result of operation at the elastic limit in 
most of the disk. More practical disks can be obtained by 
assuming <r„ values yielding greater factors of safety. The 
subject of allowable design stresses, however, requires con- 
siderable experimental investigation. 



Fig tire 5.— Stress distributions and disk profiles for illustrative designs with several assigned 
values of radial rim stress. 


SIMPLIFIED METHOD OF STRESS ANALYSIS OF GIVEN DISK 
ELASTIC STRESS 


A method for calculating the elastic-stress distribution in 
a disk of given design is presented in reference 3. Although 
this method can be directly used to check, the stress distri- 
bution in any disk designed by the procedure described 
herein, simplifications can be made that appreciably reduce 
the amount of labor involved in calculation. 

Simplification of equations. — Briefly, the method of refer- 
ence 3 consists in rewriting stress equations (1) and (2) in 
finite-difference form and solving these equations to obtain 


•ft ,n = A r , n <T t , a + B r , „ 
V t,n = A t , n (X i.a-h-Bi.n 


(14) 


where <r t>0 is the tangential stress at the first point station in 
the disk and remains unknown until the end of the calcula- 
tion, The coefficients A r n and so forth are determined from 
the equations 


A r , n — K n A r , n -iA B n At, a ~ i 

At, b=A nAf'ji- iH - A nAt,n — i 
B r ,n — KnB r ,n-\-\-L n Bt 4 n-i-{-Mn 
Bt.n=K' n B t ,n~i-{- A'nPi.n-i+A/'n „ 


(15) 


Thus, if the coefficients A„ A h B r , and B t are known at the 
(n-l)st'station, they are readily determined at the nth station 
from equations (15). The coefficients at the first station 
(r=a) are known. For a solid disk 


A r ,a — At, a — 1 
B r , a =B t , „=0 


(16) 


and for a disk with a central hole 


A r , a — B r , a — Bt , 8 — 0 ) 



Hence, successive applications of equations (15) yield the 
coefficients for subsequent substitution in equations (14) for 
actual stress determinations. The value of <T t , a is determined 
from the condition at the rim where the radial stress is known. 
Thus, from the first of equations (14) 


G t,a~ 


*r t b 




A r , i 


(18) 


where 0> i6 is the radial rim stress and A r , b and B r , t are the 
coefficients as determined in the last of the successive applica- 
tions of. equations (15). 

The bulk of the computational work in a stress analysis is 
contained in the determination of the coefficients K n , L n , and 
so forth, in equations (15). When the stations are randomly 
chosen, the equations for these coefficients are complicated. 
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A great reduction in the amount of labor is effected if the 


stations are chosen according to an assigned value of 






In addition, if Poisson’s ratio is assumed constant as n= 0.333 
the expressions assume the simplified form given in the 
following table : 


TABLE A 



050 

rn 

^ -0.200 
f* 

JK. 

0.959 t^-O-OO" 

ft* iSm-l 

o.sn^-o.ois# 2 - 

Ah ha-1 

Ln 

0.025^+0.021—- 

tin ha-I 

0.104 ^P'-rO.OH jF— 
A* -C-n-L 

K’n 

0.340 0.291 

Ah. 

0.343 0.153 ~ — 

Ah -Oh-1 

L\ 

0.009 ^‘+0.942 

Ah 

0.043 % Lt -i-0.767 

A* £ra-i 

Mn 

0.023 ~ P.-i-r0.025p») iVo*- 

- ( 0.06? -£! p„-t+0.104p„ ) r„»&P- 


0.02 iEnH'n 

0.092E.HV 

M'n 

-^0.008 P.-i-f0.009p.) r.»u'- 

-(o.02?^=- l p,-.+0.043p»^ r,«oP- 


0.977E«/r* 

0.920 EnH'n 


Alinement charts. — Evaluation of the coefficients K % , 
K'„, and so forth from the expressions given in table A is 
possible. A further reduction in labor, however, can be 
effected by the use of the alinement charts in figure 6. For 
example, in figure 6(a), joining with a straight line the points 

Eh ^ 

■^ J1 -=0.9S7 and -^=1.185 yields intersections on the suit- 
able scales for values Z*=0.199, L' x = 0.808, iT»=0.972, and 
K' n = 0.256, which check closely with calculations based on 
the expressions in table A. If the density can be considered 
constant between two stations (p*=p»_ i), the values of M x 
and M\ can also be obtained from the charts in figure 6. 

As an example, if^ =i = 1.185, r„ s pw 2 = 20.5X10 3 , and EJl r 
5.5 X10 3 , the corresponding values of and r„ a pco s are 

joined by a straight line to determine the supplementary 
point P, which is then joined to Z' a ff , ,=5.5X10 s to yield 
intersections M n =— 4.3 X10 3 and M' n =— 6.6 X10 3 . Be- 
tween any two stations where a change in density occurs, 
the expressions for M n and M\ in the table must be used. 

Illustrative example. — The necessary calculations are 
shown, in table II for the determination of the elastic stresses 
in a disk of given design, which are the disks formerly de- 
signed for the various assumed radial rim stresses. The 
coefficients K n , and so forth listed in columns 20 to 25 
are determined from values listed in table I used with the 
alinement charts of figure 6. Values beyond the range of 
the charts given in figure 6 may be computed from the 
equations given in table A. Columns 26 to 29 list the A and 
B coefficients as determined from equations (15) in which the 


values listed in equations (16) are used as the initial condi- 
tions. The value of o-, >c is then determined with the aid of 
equation (18), and the radial and tangential stresses at all 
stations are obtained from equations (14) and listed in col- 
umns 30 and 31, respectively. The equivalent stresses at 
various stations are listed in column 32, and the yield stresses 
at each station, which were the design criterions, are listed 
in column 33. A correlation of satisfactory engineering 
accuracy exists between the computed stresses and the design 
stresses shown in columns 14 and 15 of table I. 

PLASTIC STRESSES 

Modification of elastic-stress equations. — A method for 
calculating the stress distribution in a disk subjected to plastic 
flow is presented in reference 4. The form of the calculation 
is similar to that of an elastic-stress determination, except 
that the term H\ in the expression for M n and M' n (table A) 
is replaced by II' x -P' n , where P' n depends on the plastic 
strains in the radial and tangential directions. In practice, 
the plastic strains are estimated and a stress distribution is 
determined. Comparison of the plastic strains based on the 
calculated stress values with the plastic strains assumed in 
order to determine these stresses permits better subsequent 
estimates. When resultant strain values agree closely with 
assumed strain values, the strains are assumed to be correct. 

The values of P'„ to be used in the application of the sim- 
plified method are given in the following table: 


TABLE B 


Tit fji-l 
Tn 

P'n 

0.200 

0.100 A-_.-f0.125 .ir.—,-! .100Ai.»-fO.S75 

- 0.050 

0.025 Ar.n-i-0.026 At.h-i — 1.025A*,n-f-Q.974 


where A,,* is the plastic radial strain at the nth station, and 
the other A terms represent other plastic strains according 
to the notation of this report. The detailed procedure for 
estimating the plastic strains and checking with computed 
values is given in reference 4. 

Similarly for creep calculations, Q’ x terms are added to the 
P' K terms, the form of the expressions being identical to 
those for P' n in table B except that the A terms are replaced 
by creep terms S, as discussed in reference 4. 

Illustrative example. — As an illustration of practical pro- 
cedure, the disk previously designed for a radial stress of 
30,000 pounds per square inch at the rim is considered. As 
pointed out in the design calculations, any radial stress above 
0 pounds per square inch at the rim will cause the equivalent 
stress to exceed the elastic limit; hence, plastic flow must 
occur. The amount of plastic flow in this case is small and 
is confined to the single station at the rim, but the method 
involved is essentially the same for the more complicated 
cases involving plastic flow at several stations and previously 
accumulated plastic flow and creep. Examination of refer- 
ence 4 will indicate the procedure in the more complicated 
cases. 
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(a) l \ — -' -0.200. 
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Figure 6.— Alinement charts for determining constants for finltc-differeiice equations of stresses in. turbine disks. p n** p n-t. 
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The stress-strain curve for the disk material at the oper- 
ating temperature of 1070° F is given in figure 7. Point A 
is the calculated elastic equivalent stress at the rim, B the 
yield point, and C the intersection with the stress-strain 
curve of a vertical line drawn through A. In accordance 
with the method of reference 4, the strain CD is the first 
approximation of the plastic strain. The value is listed in 
columns 34 and 35 of table III for station b. Columns 36 
and 37 list the approximated radial and tangential plastic 
strains obtained from equations (20) of reference 4, The 
value of P' n as determinecLfrom. table B, is listed in column 
38. Column 39 corresponds to column 8 of table I used in 
the elastic-stress calculations; the only difference is that 
fiT'-P' has replaced H' n . Prom this point-, the calculations 
are identical to elastic-stress calculations and are shown in 
columns 24, 25, and 28 to 32 (table III), which are similar 
to the corresponding columns of table II. Further approxi- 
mations, made in accordance with the method of reference 
4, are also shown in table III to obtain greater accuracy of 
approximation. 


COMMENT ON SIGNIFICANCE OF THERMAL STRESSES 

The design method described assigns equal importance 
to the thermal and centrifugal stresses. The true relative 
importance of these two types of stress is, however, as yet 
to be evaluated. Whereas the centrifugal stresses are 
inherently necessary to balance the centrifugal forces of 
rotation— so that yielding due to overstressing in one region 
shifts the load and thereby increases the stress in another 
region— the thermal stresses constitute an internally balanced 
system. Tension in one region balances compression in 
others; hence, plastic flow in one region reduces the thermal 
stresses in another region. When considering bursting 
strength, for example, centrifugal stresses are probably of 
greater” importance than thermal stresses. In the case of 
rim cracking (reference 8), which is in part due to cyclic 
plastic, flow at the rim, the thermal stresses are, however, 
fully as important as the centrifugal stresses. Likewise, 
when disk deformations are considered in designing for close 
mechanical tolerances, thermal stresses must receive atten- 
tion equal to that given to centrifugal stresses. 

The relative importance of centrifugal and thermal 
stresses depends, to a great extent, on the ductility of the 
disk material. As yet, experimental data are very meager 
on the significance of ductility. In any case, however, 
because the distinction between the two types of stress can 
only be realized after the occurrence of plastic flow and 
because the behavior of a material in subsequent loading 
depends on its previous stress-strain history, thermal stresses 
should receive detailed attention early in the design in order 
to minimize any detrimental effects that may result from 
their omission. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 4 , 1949. 
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In addition to the symbols used herein, the symbols 
of references 3 and 4 of which application is made in this 
report are defined. The various coefficients are given in 
their general form for random spacing of the point stations. 
From these coefficients, the simplified coefficients can be 


obtained by direct substitution of 
p=constant. 


r„— r K _ t 


= constant and 


base area of airfoil section of blade, sq in. 

E elastic modulus of disk material, lb/sq in. 
h axial thickness of disk, in. 

N number of blades supported by disk 
r radial distance, in. 

T temperature, °F 

a. coefficient of thermal expansion between actual 
temperature and temperature at zero thermal stress, 
m./(in.)(°F) 

A plastic increment of strain, in./in. 

A r plastic increment of strain in radial direction, im/in. 

A t plastic increment of strain in tangential direction, 
in./in. 

A T temperature increment above temperature of zero 

thermal stress, °F 

S r creep increment in radial direction, in./in. 

S t creep increment in tangential direction, in./in. 

e p plastic strain corresponding to stress ay in tensile 
specimen, in./in. 

/ 1 Poisson’s ratio 

p mass density of disk material, (lb) (sec*) /in. 4 
o-j i stress at base of airfoil section of blade, lb/sq in. 

ay equivalent stress, lb/sq in. 

ay radial stress, lb/sq in. 

ay tangential stress, lb/sq in. 

ay proportional elastic limit, lb/sq in. 

a angular velocity, radians/sec 

The following supplementary subscripts are used for de- 
noting values of the proceeding symbols in connection with 
the finite-difference solution: 


nth point station 
(7i-l)st point station 

station at smallest disk radius considered (For disk 
with a central hole, this station is taken at the 
radius of the central hole; for a solid disk, this 
station is taken at a radius approximately 10 
percent of the rim radius, which is larger than that 
(5 percent) recommended in reference 3. The 
value of 10 percent yields results of satisfactory 
accuracy and permits the use of fewer stations 


when constant value of 



is used.) 
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( M n and M' n are defined in the elastic case for 

P'n=Q'«= 0 ) 


For ^-^ =0.050, 



Z»=W*- (0.965er*. ,_t- 0.298a r . -EJI' 

■P'n-l 

For r - n —— n ~ 1 = 0.200, 
r« 




(0.833<r t , 0. 166tr r , »_i) — E n H' 
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TABLE I.— CALCULATION OF DISK PROFILES FOR ILLUSTRATIVE DESIGN 



I 

2 

3 

4 

5 

6 

7 

8 

9 

10 

It 

Tit - Tn—I 

r™ 

r«. 

AT. 

a. 

ct* 


H'„ 

EJT,. 

r »VA 

a. 

En-l 


T* 


r.XCO 




OX® 

(63-®.-! 

(43 X (7) 

(2JIX1I50 

(4) + (4).-i 

a 

0.200 

0.100 

0.S4S 

370 

28.6X10* 

8. 88X10-* 

3.288X10-* 





2 

.200 

-132 

L 056 

375 

28.6 

8.90 

3.33S 

0.049XKH 

1. 40X10* 

J. 28X10* 

1.000 

3 

.200 

.165 

L 320 

3S0 

28.5 

8 80 

3.382 

.044 

1.25 

2.00 

.997 

4 

.200 

.206 

4648 

385 

28.5 

8. 91 

3.430 

.015 

1.37 

3.12 

1.000 

5 

.200 

.258 

2. 004 

395 

28.4 

8.93 

3.527 

.097 

2.75 

4.90 

.996 

6 

.200 

.322 

2. 576 

410 

28.3 

8.95 

3.670 

.143 

4.05 

7.63 

.996 

7 

.200 

. 103 

3. 224 

435 

28.2 

9.00 

3.915 

.245 

6.91 

11.95 

.936 

8 

.200 

.504 

4.032 

470 

28.0 

9.05 

4254 

.339 

9.49 

18.70 

.993 

9 

.200 

.630 

5.010 

560 

27.4 

9.20 

5. 152 

.898 

24.61 

29. 19 

.979 

10 

.050 

.663 

5.304 

580 

27.3 

9,23 

5.353 

.201 

5.49 

32.35 

.936 

11 

.050 

.698 

5.5S4 

610 

27.0 

9.28 

5.661 

.308 

8. 32 

35.86 

.989 

12 

.050 

.735 

5.8S0 

615 

26. S 

9.33 

6.013 

.357 

9.57 

39.76 

.393 

13 

.050 

.773 

6.184 

685 

26.5 

9.40 

6. 439 

.421 

11. 16 

43.93 

.989 

14 

.050 

.814 

6.512 

730 

26.1 

9.48 

6.920 

.4S1 

12.55 

48.77 

.985 

15 

.050 

.857 

6.856 

785 

25.6 

9 56 

7.505 

.585 

14.98 

5405 

.98J 

16 

.050 

.902 

7.216 

S45 

248 

9 65 

8.163 

.638 

16.32 

59.88 

.369 

17 

.050 

.950 

7.600 

910 

23.5 

9.77 

8.S91 

.728 

17. II 

66.12 

.943 

b 

.050 

L 000 

8-000 

1000 

21.4 

9.92 

9.920 

L 029 

22.02 

73.60 

.911 



II 

12 

13 

14 

15 

16 

17 

18 

19 

ft 

0.833X(153.-i 

-0.166X(14),.-i 

z». 

(11)X(10)— (83 

OV.K 

(Tr,8 

<rt,n 

8.00X(143»-i-f- 
(153 .. 1 — 0.64XC9) 

IO.OOX(14) 
— (15)4(9) 

A»-i 

A.'’ 

(17) -(IS) 

A„ 

(IS)«iX(19).« 

0>.*~ 

18,000 

C T .* = 
20,000 

<r t . 

30,000 

Cr. i~ 

40,000 

a 



76.5X10* 

70. 5X10* 

78.5X10* 




1. 86 

1.75 

1.50 

1.38 

2 

5L0X10* 

49.6X10* 

76.3 

76.0 

75.0 

688X10* 

686X10* 

0.997 

1.87 

1.76 

450 

1.36 

3 

49.9 

48.5 

76.2 

76.0 

740 

682 

6SS 

4009 

1.85 

1.74 

1.49 

1.35 

4 

49.0 

47.6 

76.1 

77.0 

73.0 

680 

700 

L029 

1.S0 

4 69 

445 

431 

5 

48.0 

45.1 

76.0 

78.0 

72.0 

6S6 

713 

1.039 

473 

1.63 

1.40 

428 

6 

47.0 

42.8 

75-9 

79.0 

70.5 

691 

727 

1.052 

464 

455 

433 

420 

7 

45.6 

38.5 

75.5 

80.5 

66-5 

695 

750 

4079 

1. 52 

444 

423 

1. 11 

S 

42.0 

32.2 

75.0 

83.0 

62.0 

699 

787 

4126 

1.35 

428 

LOS 

.99 

9 

37.9 

12.5 

73.5 

840 

46.0 

707 

823 

L 164 

416 

410 

.94 

.85 


0-965X(153.-i 





38.00X(14),-i4- 

40.00X(14) 







-0.298X(14)^i 





(153.-i-0.903X (9) 

-(15)4-19) 






10 

19.4X10* 

13.8 

73.0 

83.5 

42.5 

! 3209X10* 

3330X10* 

1.038 

1.12 

1.05 

.91 

.82 

11 

16.1 

7.6 

72.5 

83.0 

36.5 

3183 

3319 

4043 

L 07 

1-02 

-87 

.79 

12 

10.5 

.8 

7L5 

82.0 

29.5 

3155 

3290 

1.043 

1-03 

.98 

.83 

.76 

13 

40 

-7.2 

70.5 

78.5 

21.0 

3106 

3163 

1.018 

1. 01 

.96 

.82 

.75 

14 

-3.1 

-15.6 

69.5 

745 

11.5 

2960 

3017 

1. 019 

.99 

.94 

.80 

.74 

15 

-11.1 

-25.9 

68.0 

67.0 

-1.0 

2794 

2735 

.979 

L 01 

.96 

.82 

.76 

16 

-20.9 

-36.6 

65. S 

55.5 

-16.0 

2491 

2336 

.938 

4 08 

1.02 

.87 

.81 

17 

-32.3 

-47.7 

63.0 

43.0 

-31-0 

2071 

1817 

.877 

1.23 

1. 16 

.99 

.92 

b 

-42.7 

-60.9 

60.0 

0 

-60-0 

1537 

134 

.087 





b 

—42.7 

-60.9 

63.5 

18.0 

-53.0 

1537 

847 

.551 

2.23 




b 

-42.7 

-60.9 

64.0 

20.0 

-52.0 

1537 

926 

.602 


L93 



b 

-42.7 

-60.9 

67.7 

30.0 

-48.0 

1537 

1322 

.860 



L 15 


b 

-42.7 

-60.9 

72.9 

40.0 

-44.5 

1537 

1718 

1.118 


— 


.82 
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TABLE II— CALCULATION OF ELASTIC STRESSES FOR ILLUSTRATIVE DESIGN 


From table I, 
(18) and (10) 


Coefficients determined from figure 6 
or from table A 



-a4xios 

-1. 4X10* 

-.5 

-1.3 

-.7 

-1.5 

-1.1 

-2.9 

-1.7 

-4.3 

-2.7 

-7.2 

-4,2 

-10.1 

-7.6 

-24.8 

-1.7 

-5.9 

-2.0 

-8.8 

-2.2 

-10.0 

—2.4 

-11.7 

— 2.7 

-13.1 

— 2. 9 J 

-15.5 

-3. 2 - 

-16.9 

-3.4 

-17.8 



, fMvf28 , , C20)X(28)„-r (22>X(2S).-, -76, 100 

iSSSr*, Swran , +C23)XC29) n -< ±± 

<ri..X(26) <ri,.X(27) 

+(28) +(29) 


74.6 

78.6 
73.3 

71.0 

70.0 

06.6 
62.5 
45.8 
42.7 


0 

76. 1X10* 

-1. 4X10* 

76.4 

-2.5 

75.3 

-3.7 

76.3 

-6.3 

77.3 

-10.1 

78.3 

-16.7 

79.9 

-25.9 

82.4 

-49.7 

83.5 

-54.8 

83.1 

-62.5 

82.7 

-71.5 

81.7 

-81.4 

78.2 

-92.2 

74.2 

-104.0 

60.8 

-114.7 

56.4 

-122.4 

42.9 



TABLE III.— CALCULATION OF PLASTIC STRESSES FOR ILLUSTRATIVE DESIGN 

[»,.!«» 30,060] 
























































